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Abstract 

It is found that the evolutions of density perturbations on the super-Hubble 
scales are unstable in the model with dark-sector interaction Q proportional to 
the energy density of cold dark matter (CDM) p m and constant equation of state 
parameter of dark energy Wd- In this paper, to avoid the instabilities, we suggest 
a new covariant model for the energy-momentum transfer between DE and CDM. 
Then we show that the the large-scale instabilities of curvature perturbations can 
be avoided in our model in the universe filled only by DE and CDM. Furthermore, 
by including the additional components of radiation and baryons, we calculate the 
dominant non-adiabatic modes in the radiation era and find that the modes grow 
in the power law with exponent at the order of unit. 

PACS: 95.36. +x, 98.70.Vc, 98.80.Cq 
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1 Introduction 

We are convinced by the increasing observations [H EJ [3] that the present universe is 
dominated by the two components: dark energy (DE) [HE] and cold dark matter (CDM). 
Currently, the two components are only indirected detected via the total gravitational 
effects. And then this produces an degeneracy: the two dark components might interact 
mutually without violating the observational constraints [6l [7] . And furthermore it is 
found that an appropriate interaction can help to alleviate the coincidence problem [8l 
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namely why DE and CDM are comparable in size exactly today [5]. Different interacting 
models of dark energy have been investigated intensively [TOj ITT] . 

Generally, a model with interaction between DE and CDM is described in the back- 
ground by the two continuity equations 



where Q denotes the phenomenological interaction term between DE and CDM; p m and 
Pd are the energy densities of CDM and DE respectively; Wd = Pd/ Pd is the equation 
of state parameter of DE; pd is the pressure density of DE; H = a/ a is the Hubble 
parameter; a(t) is the scale factor in the Friedmann-Robertson- Walker (FRW) metric; 
a dot denotes the derivative with respect to the cosmic time t. In the note we do not 
allow the phantom case Wd < — 1. Owing to the lack of the knowledge of micro-origin of 
the interaction, usually the interaction term is parameterized in a simple form as [8] 



where a and f3 are positive constants. The interaction term Q would influence not only 
the background dynamics of the universe, but also the growth of the perturbations of 
the cosmological fluids [T5J EH [33] • 

Recently, in Ref.|15j. by modeling DE as a fluid with constant Wd > — 1, the au- 
thors investigated the evolution of the linear density perturbations and found that the 
combination of constant Wd and the simple interaction form Q given in Eq.Q leads to 
an instability: the curvature perturbation on the super-Hubble scales blows up in the 
early universe [T5]. The explicit models investigated in Ref.[T5] included the two cases of 
(3 = and a = (3 in Eq.([n]). Further more, in [17], it is concluded that the perturbations 
in the dark energy become unstable for any model with constant Wd > — 1 and non-zero 
a, no matter how small the parameter a is made. In [161 H3 HB], t ne case of a = was 
surveyed and it was found that the instability can be avoided if (3 is made small enough. 
In [TU], by modeling DE as a quintessence field, the author found that the instability 
can also be avoided even for the interaction Q proportional to p m . 

Then it seems that a model with constant Wd and non-zero a in Eq. ([3]) would be ruled 
out as a viable interacting model. To cure the curvature perturbations, the authors in 
Ref.[20J suggested a covariant model for energy-momentum transfer between DE and 
CDM, and showed that the instabilities of the curvature perturbation on the large scales 
can be avoided in the covariant model with constant Wd and Q = 3aHp m . In Ref.[20j, 
the authors defined the effective EMTs of CDM and DE respectively as 



Pm + 3Hp m = Q, 
p d + 3H(1 + w d )p d = -Q, 



(1) 
(2) 



Q = 3H(ap m + (3p d ), 



(3) 



1 em rm ll m u m ' t , m\ a m u "m ' i) ) 



(4) 
(5) 
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where p^ = —ap m , p c d = pd + ap m , and and are the four velocities of CDM and 
DE respectively. The two effective EMTs are taken to be conserved [20] 

TZ-,u = 1ST.,, = 0, (6) 
and the Einstein equation are assumed to be [20J 

CT = BnG(TZ + I%). (7) 

Based on Eqs.© and (J7|), the evolution of the curvature perturbations were surveyed 
in [20] and no instabilities were found. The model defined in Eqs.(|6]) and (JTj) has a 
distinguishable feature that the total EMT of DE and CDM is not conserved at the 
perturbative level (Tgf + T^ v \ v ^ 0, while, in other covariant models, the total EMT 
is always taken to be conserved in order to match the Einstein tensor. Of course, in 
[20] , the non-conserved total EMT does not cause problems since it is the total effective 
EMT that appears on the right-hand side of the Einstein equation and is guaranteed to 
be conserved by Eq.fjSJ). 

But it is still puzzling and discomforting that the total EMT of DE and CDM is 
not conserved. We try to solve the puzzlement by suggesting a new covariant model in 
this paper. We consider DE as a fluid with constant Wd that is coupled to CDM via a 
covariant energy-momentum transfer, from which Q = 3aHp m can be reduced at the 
background level. We show that the total EMT of CDM and DE is conserved and the 
stabilities of the curvature perturbation on the large scales can be avoided in our new 
covariant model. 

This paper is organized as follows. In SecEl we display our new covariant model for 
the energy-momentum transfer in the dark sector. In Secj3], by assuming the universe 
filled only by DE and CDM, we survey the evolution of the density perturbations and 
show that instabilities on the large scales can be avoided. In SecJU by considering the 
effects of the radiation (photons and neutrinos) and baryons, we investigate the dominant 
non-adiabatic mode in the radiation era and show that no non-adiabatic modes blow up. 
Finally, conclusions and discussions are given. 



2 New Covariant Model for Interaction in Dark Sec- 
tor 

Usually, the covariant form for energy-momentum transfer is taken to be [151 ETJ 

V„Tf = Q% (8) 

where A — m, d to denote CDM and DE respectively, and the condition Y2a Qa = is 
imposed in order for the total energy-momentum tensor to be conserved. By comparing 
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Eqs.flHD an d ®, we can get the covariant interacting terms of the model in Eq.® as [20] 

Q^ m = a[p m (u^ m + 9 nU (9) 
Q» = -a[p m (u»u» d + 9n];v (10) 

Clearly, at the perturbation level, Q^ n ^ —Qd and so the total EMT of DE and CDM is 
not conserved. 

Yet, motivated by Eqs.flH]) and f JTOj) . we find that a new covariant model can be 
constructed as 

T£, v = (P m «j;, = ~Qd, (11) 

n v , v = [ Pd <u» d + p*w + <r)]v = (12) 

where Q d is given in Eq. (fTU|) . Obviously, in this model the total EMT of DE and CDM 
is conserved. Indeed, this new model is just a usual type of Eq.flS]) with the choice of 
the covariant energy-momentum transfer as 

Our model is similar to the case in [22| |23| 124] where 

Qm = -Qi = mv^vh 

in the sense that the interaction is determined by the energy density of CDM and the 
four velocity of DE. It can be easily checked that Eqs.((T]) and (J2J) with Q = 3aHp m 
can be deduced from Eqs.f TTT]) and ( TI2"]) at the background level, respectively. Notably, 
this model has a similar feature to the one in Ref.JSU] that the global quantity H in 
Q = 3aHp m is explained to be a local quantity u v d . v . 
Now together with the Einstein equation 

G^ v = 8ttG(T^ + T£ v ), (13) 

we can study the evolution of the curvature perturbation on the large scales to check 
whether the large-scale instabilities could be avoided. 

3 Evolution of Density Perturbations 

In this section, we apply the new model suggested in the last section to survey the 
evolution of the density perturbations by modeling DE as a fluid with constant Wd- 
For simplicity, we consider a flat universe filled only by DE and CDM. We choose the 
conformal Newtonian gauge and then the perturbed FRW metric in the conformal is 
given by 

ds 2 = a 2 (r)[-(l + 20)rfr 2 + (l-2^)c/x 2 ], (14) 
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where 4> an d ip denote the scalar perturbations. The Friedmann equation reads 

n 2 =\-\ =—( Pm + Pd ) (15) 



Hereafter, primes denote the derivatives respect to r. With Q = 3aHp m and constant 
w d , we can obtain the background evolutions of p m and pd from Eqs.flT]) and (J2J) as 

Pm=Pm a- 3{1 - a \ (16) 

Pd = Pdoa- 3 ^ 1+Wd) + (^—) Pm »a-\a-^ - a 3 "). (17) 
V a + w d ' 

In this paper, we use the subscript to denote the present value of the corresponding 
parameter and a = 1. 

3.1 Evolving equations of Perturbations 

When the perturbed metric in Eq. (|14p is considered, the four velocities of CDM and DE 
are 

u£ = a~ x fl - (j), diVrn) , u^ d = a _1 (l - (j),diV d \ (18) 

where v m and v d are the peculiar velocity potentials of CDM and DE respectively. Usu- 
ally, we define the volume expansion rates of CDM and DE (in Fourier space) respectively 
as 

e m = - k 2 v m , e d = -k 2 Vd . (19) 

We use Sp m , 5p d and 5p d to denote the first-order perturbations of the corresponding 
parameters, and introduce two dimensionless first-order parameters as 

r S Pm x Spd {on , 

d m = , Od = • (20) 

Pm Pd 

The total curvature perturbation on the constant-p (p = p m + p d ) surface is defined as 

( = -^-U 5 -^±^. (21) 

P'm+P'd 

Here, following the analysis in Ref.[15j, we take 

Spd = Sp d + (1 - w d )[3H(l + w d + a— )p d ]%. (22) 

Pd k 

Then from Eq. fTTB"]) . we get two equations as 

S' m - 3(1 - a)V>' + 9 m - a9 d = 0, (23) 
6' m + H{\ + 3a)6 m - k 2 (l -a)<f> = a [6' d + H(l + 3a)6 d - k 2 6 m ] . (24) 
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And from Eg. ([IT]) , we get other two equations as 

5 > + 3W _ Wd _ a P^) 5d + 9^2 (1 _ +Wd + a ElIL)^ 

Pd Pd k 2 

+ 3aU—6 m + (1 + w d + a— )(0 d - 3^') = 0, 

a(l — a) 2 ^ 

Q' d - 2U6 d - m- &-0 d - k 2 <\> 

d l + w d + a^ 

Pd 

u2 Pi -S m - k 2 °-* ^ = 0. 



(25) 



(26) 



1 + w d + m 1 + w d + 

Pd a pd 



In the conformal Newtonian gauge, the first-order Einstein equations (TT3|) gives us 

m 

3W + fcV + 3^ 2 = -4nGa 2 (5 mPm + 5 d p d ), (27) 

fcV + fc 2 ft<£ = 4vrGa 2 [p m ^ m + (1 + w d )pA] (28) 

V>" + W(2^' + 0') + (2— - U 2 )<\> +%.($-<!>) = A-KGa 2 5p d (29) 

a 3 

V> - = 0. (30) 

Only two of the above equations are independent. Choosing any two of them and using 
Eqs. (|23p - (|26p . we can solve these evolving equations numerically if the initial conditions 
are given. 



3.2 Adiabatic Initial Conditions 

In the early universe, a < 1, Eqs. flTB"]) and (fT7j) indicate 



Pm w d + a 
Pd a 



and then, from Eq. f!15j) . we have 



= T^3a ' r= (l-3^oV^ a ' ^ 

where 

= 87rGp m0 
m0= 3# 2 ' 

Here we adopt the adiabatic initial conditions to study the evolution of the density 
perturbations on the super-Hubble scales ik <C aH). To the lowest order in kr, we can 
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Fig. 1: log 10 |C| versus log 10 a in the interacting model for fixed Q m o = 0.3, h = 0.67, 
k = 1.5 x 10~ 4 Mpc _1 , a = 10~ 3 and different Wd- 

set the adiabatic conditions as 

= if) = = Const., (33) 
5 m = 5 d = -2A+, (34) 

m = 9 d = 1-3a fc 2 rA . (35) 
3(1 — a) 

3.3 Evolution of Curvature Perturbation on Large Scales 

Now using the initial conditions in the last subsection, we can obtain the evolution of 
the scalar perturbations by solving Eqs.( l23|) -( l27j) and (130]) numerically, and then get the 
evolution of the curvature perturbation ( defined in Eq. (T2Tj) . We display the results in 
Figfj] and FigfSJ Here we have taken = 10~ 25 and a = 1, and fixed the values of 
some parameters: fi m o = 0.3, k — 1.5 x 10 _4 Mpc _1 , Hq = lOO/i km sec _1 Mpc _1 and 
h = 0.67. The evolving curves in Figfj] and Figj2] manifest the regular growth in power 
law and no instabilities occur. 



4 Dominant Non-adiabatic Modes 

In the last section, by solving the evolving equations numerically with the adiabatic 
initial conditions, we show that instabilities of the curvature perturbation on the large 
scales are removed in our new covariant model. But the conclusion is obtained in the 
universe filled only by CDM and DE. In this section, we consider the effects of baryons, 
photons and neutrinos by calculating the dominant non-adiabatic modes deep in the 
radiation era. If the dominant non-adiabatic modes do not grow rapidly, we believe that 



7 




Fig. 2: log 10 |C| versus log 10 a in the interacting model for fixed fl m0 = 0.3, h = 0.67, 
k = 1.5 x 10~ 4 Mpc -1 , w d = -0.94 and different a. 



the instabilities can be avoided even when radiations and baryons are involved in the 
universe. 

The A-fluid EMT including perturbations is taken to be 

= (pa + Va)u\u a + PAS T + <\ (36) 

where u A is the four velocity 



u A = a 1 [ 1 



We have allowed an anisotropic shear perturbation 7r^, and A = m,d, b, 7, v to denote 
the corresponding parameters of CDM, DE, baryons, photons and neutrinos. Following 
Ref.pS], we take = iff = 7if = tt^ = and 

71^ = 0, n'J = a- 2 (d i d j -^)Tr v . (37) 

For DE and CDM, the continuity equations fill I) and 0121) still hold. 
Early in the radiation era, the Friedmann equation reads 

2 n 8tiG 8ttG _a 

Ha = —^-{P-y + Pu) = —^-Proa ■ (38) 

Here we use the subscript r to denote the corresponding parameter of radiation which 
consists of photons and neutrinos. Then we have 



Uoflor, h = t-\ Q r0 = ^M. (39) 



In the radiation era, the perturbed Einstein equations give us that 

3r - y + k 2 ifj + 3r" 2 = -4vrGa 2 (8 Pm + 5p d + ^ <5p A ) , (40) 

/c 2 y + r~V) = 4vrGa 2 p m # m + (1 + w rf )p d ^ + ^(p A + Pa)#a] , (41) 

A 

k 2 f \ 

i>" + 2r~y + r _ V " T~ 2 (j) + T {ip - (j>) = ^Ga 2 [5p d + V 8p A ) , (42) 

i/>-(f> = 8nGir u , (43) 

where A runs over b, 7, and za 

For CDM and DE, the perturbed continuity equations are given by Eqs. 
For baryons, the perturbed continuity equations (in Fourier space) are 



and for photons [15] 



5' b = -9 b + 3^', (44) 
9' b = -W b + /c 2 0, (45) 



5; = ~0 7 + 4^, (46) 
0; = ^ 2 5 7 + /c 2 0, (47) 



and for neutrinos 



8' v = -\e v + hl>\ (48) 
^ = ifc 2 5, + k 2 <P - k 2 a u , (49) 

< = 7^' (50) 
15 

where 0^ = —k 2 VA for A = 6, 7, z/, and cr^ = 2/c 2 7r ; ,/[3a 2 (p i , + p^)] 

In order to find the dominant non-adiabatic modes, we assume a leading-order power 
law for perturbations [15] 



V = A^r) n *, = ^(Arr)"*, 5 a = B A {kr) n \ 6 A = C A (kr) s \ a v = D v (kr) n " . (51) 

Here the subscript A = m,d, b, 7, and z/ denotes the corresponding parameter of CDM, 
DE, baryons, photons and neutrinos, respectively. To the leading order in kr, the equa- 
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tions (|23l) - (l26j) and (I40l - (l50l) may be solved, in terms of if): 



6 = Jib, J = 1 - w 16jR " , (52) 

V 5(?v + 2)(r^ + l) + 8iV 1 ; 

5^ = 5 U = Aif), h = 3-0, (53) 

7-1-1 4 rf) 

e J = e u = e b = ——k 2 r^ } o u = --^- (54) 

n^, + 1 15n^ + 2 

5 m = 3(1 - a)^, (55) 
2ft (1 ~ 3o)/2 i/> 

#d = —ZT7Z — Ui ' ttt rfc 2 r5 d , 6> m = a6> d , (57) 

9(1 - iu d )(l - a) 

where R u = p^/ (p 7 + p^) and 

-3 Wd ± y/9wj + I2w d - 20 
rty = , (58) 

From Eq.f l56|) . we should require 

Re[n^] > 1 + 3a. 

in order for the modes to be well behaved as kr — > 0. For wa ~ — 1, this leads to 

3 / 1 

--Wd > 1 + 3a => a < -. 

Correspondingly, the total curvature perturbation ( is defined as 

I^aPa 

where A runs over m, d, 7, v and 6. Then £ can be expressed in terms of if) as 

C = ~K + J + 2)V- (59) 
For ~ — 1, is a complex number and 

t. r 1 3 

So the dominant non-adiabatic modes grows in a power law with exponent at the order 
of unit and no instabilities are present. 
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5 Conclusions and Discussions 



It was found in [15] that the evolution of the density perturbations is unstable in the 
model with constant Wd and Q proportional to p m . To cure it, a covariant model for 
energy-momentum transfer in dark sector was suggested in j2U]. Yet, the model in [2U] 
seems to be "anomalous" , since it has a puzzling and discomforting feature that the total 
EMT of DE and CDM is not conserved. In order to remove the "anomaly" , in this paper, 
we have suggested a new covariant model in which the total EMT of DE and CDM is 
conserved. By choosing Q = 3aHp m and using the new model, we survey the evolution 
of density perturbations in a universe filled by DE and CDM, and show that the large- 
scale instabilities can be avoided. Furthermore, we calculate the dominant non-adiabatic 
modes in the radiation era by including the other components of radiation and baryons 
in the universe, and find that the dominant non-adiabatic modes grow in the power law 
with exponent at the order of unit. This makes us believe that the instabilities can also 
be avoided even in a universe filled by radiation, matter and DE. 

Qualitatively, it is easy to understand why the instabilities can be avoided in our 
new model. In the model in [15] . it was shown that the coupling term Q leading the 
instability-driving term in the perturbed continuity equation for 9d during the radiation 
era 

ff d „ — ^ — —TiOd = -^W, (60) 

If Wd is close to —1, this term becomes very large and causes the rapid growth of 9 d 
during the radiation era. Nevertheless, in the model in J2U], it was shown that the 
corresponding driving term in the continuity equation for 9d during the radiation era 
becomes 

1 + w d + a{p m /pd) Pd I- a 

Clearly, this term does not cause instabilities even for Wd ~ — 1. So it is believed that 
the instabilities can be avoided in the model in [20]. In fact, the continuity equation for 
9d in our new model, Eq.( l26|) . are same to the one in [20]. Then the analysis based on 
Eq.( 1oTj) also apply here. So we believe that the instabilities can be avoided in our new 
model, as in the model in 
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